SIMULATION AND APPROXIMATION OF LEVY-DRIVEN 
STOCHASTIC DIFFERENTIAL EQUATIONS 



Abstract. We consider the problem of the simulation of Levy-driven stochas- 
tic differential equations. It is generally impossible to simulate the increments 
of a Levy-process. Thus in addition to an Euler scheme, we have to simulate 
approximately these increments. We use a method in which the large jumps 
are simulated exactly, while the small jumps are approximated by Gaussian 
variables. Using some recent results of Rio 1121 about the central limit the- 
orem, in the spirit of the famous paper by Komlos-Major-Tsunady | 10 |. we 
derive an estimate for the strong error of this numerical scheme. This error 
remains reasonnable when the Levy measure is very singular near 0, which is 
not the case when neglecting the small jumps. 

In the same spirit, we study the problem of the approximation of a Levy-driven 
S.D.E. by a Brownian S.D.E. when the Levy process has no large jumps. 
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1. Introduction 

Let (Z t )t>o be a one-dimensional square integrable Levy process. Then for some 
o£l, some b £ R+, and some measure v on M* satisfying L z 2 v{dz) < oo, 



where (B t )t>o is a standard Brownian motion, independent of a Poisson measure 
N(ds, dz) on [0, oo) x R„ with intensity measure dsv(dz), and where N is its com- 
pensated Poisson measure, see Jacod-Shiryaev [9]. 

We consider, for some x £ E and some function a : R ^ K, the S.D.E. 



Using some classical results (see e.g. Ikeda-Watanabe [4]), there is strong existence 
and uniqueness for ([2]) as soon as a is Lipschitz continuous: for any given couple 
(B,N), there exists an unique cadlag adapted solution (X t )t>o to ([2|). By adapted, 
we mean adapted to the filtration (Ft)t>o generated by (B, N). 

We consider two related problems in this paper. The first one, exposed in the 
next section, deals with the numerical approximation of the solution (Xt)t>o- The 
second one concerns the approximation of (Xt)t>o by the solution to a Brownian 
S.D.E., when Z has only very small jumps, and is presented in Section[3] 
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Our results are based on a recent work of Rio [12] that concerns the rate of conver- 
gence in the central limit theorem, when using the quadratic Wasserstein distance. 
This result, and its application to Levy processes, is exposed in Section 2J 

2. Numerical simulation 

The first goal of this paper is to study a numerical scheme to solve {2} . The first idea 
is to perform an Euler scheme (Xy n ) n > with time-step 1/n, see Jacod [5], Jacod- 
Protter [8], Protter-Talay [11] for rates of convergence. However, this is generally 
not a good scheme in practise, unless one knows how to simulate the increments of 
the underlying Levy process, which is the case e.g. when Z is a stable process. 

We assume here that the Levy measure v is known explicitely: one can thus simulate 
random variables with law v((1z)1a(z)/v(A), for any A such that v{A) < oo. 

The first idea is to approximate the increments of Z by A™' 6 = Zf, — ^fi—iy n > 
where Z\ is the same Levy process as Z without its (compensated) jumps smaller 
than e. However, Asmussen-Rosinski [Ij have shown that for a Levy process with 
many small jumps, it is more convenient to approximate small jumps by some 
Gaussian variables than to neglect them. We thus introduce A"' e = A"' e + £/™' e , 
where U™' e is Gaussian with same mean and variance as the neglected jumps. The 
arguments of [1] concern only Levy processes, and it does not seem so easy to apply 
such an idea to the simulation of SDEs. 

Let us write (X^y n ) t >o (resp. (^^i/„)i>o) f° r the Euler scheme using the ap- 
proximate increments (A™' e )i>i (resp. (A™' e )i>i). They of course have a similar 
computational cost. 

Jacod-Kurtz-Meleard-Protter [7] have computed systematically the weak error for 
the approximate Euler scheme. In particular, they prove some very fine estimates 
of E[j(lK, )] — E[<?(Xt)] for g smooth enough. The obtained rate of convergence 
is very satisfying. 

Assume now that the goal is to approximate some functional of the path of the 
solution (e.g. supj r i Then we have to estimate the error between the laws 

of the paths of the processes (not only between the laws of the time marginals). 
A common way to perform such an analysis is to introduce a suitable coupling 
between the numerical scheme (X^y n )t>o and the true solution (X t )t>o, and to 
estimate the (discretized) strong error E[sup[ T ] \X^, — X[ nt ]/ n | 2 ]. 
We refer to Jacod- Jakubowski-Memin [B] for the rate of convergence of the dis- 
cretized process (X[ nt y n ) t >o to the whole process (X t ) t >o- 

Rubenthaler [13] has studied the strong error when neglecting small jumps. He 
obtains roughly E[sup [0 T] \X^ ]/n - X [nt]/n \ 2 } ~ C T {n~ l +Jj z |< e z 2 v(dz)) (if 0). 
For v very singular near 0, the obtained precision is very low. 

Our aim here is to study the strong error when using X"^y n . We will see that the 
precision is much higher (see Subsection 12.41 below). 

The main difficulty is to find a suitable coupling between the true increments (Z i j n — 
Z(i-i)/ n )i>\ and the approximate increments (A™' e )i>i: clearly, one considers Z, 
then one erases its jumps smaller than e, but how to build the additionnal Gaussian 
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variable in such a way that it is a.s. close to the erased jumps? We will use a 
recent result of Rio |12j , which gives some very precise rate of convergence for the 
standard central limit theorem in Wasserstein distance, in the spirit of Komlos- 
Major-Tsunady [TO] . 

2.1. Notation. We introduce, for e € (0, 1), k £ N, 

(3) F e (y) = f v{dz\ m k {v) = f \z\ k v{dz), 



m fc »= / \z\ k is(dz), 8 e (v)- '" L,{ " ] 



\z\<e m2,e(v) 

Observe that we always have S^v) < e 2 and F e (y) < fT 2 mi(y). 

For n £ N and t > 0, we set p n (t) = [nt]/n, where [x] is the integer part of x. 

2.2. Numerical scheme. Let n E N and e £ (0, 1) be fixed. We introduce an i.i.d. 
sequence (A"' e )i>i of random variables, with 

N„, € 

(4) A^ e = a n , e + b n>e G+J2Yi, 

i 

where a„ i£ = (a — J| z | >£ zv(dz))/n, where b\ e = (b 2 + m2. e {v))/n, where G is 
Gaussian with mean and variance 1, where N n ^ is Poisson distributed with mean 
F e (v)/n, and where Y^Kf,... are i.i.d. with law ^(<iz)l| 2 | >e /i 7 ' £ (^). All these 
random variables are assumed to be independent. Then we introduce the scheme 

(5) X^ = x, ^ 1)/n =^+a(l-)A- «>0). 
Observe that 

• the cost of simulation of A™' e is of order 1 + E[JV ni€ ] = 1 + F e {y)/n, whence that 
of ( X p£t))telo,T\ is of order Tn(l + F e {v)/n) = T(n + F e {v)), as in [J3]; 

• A™^ has the same law as Zf i+1 y n — Zf, n + U n , e , where U n ^ is Gaussian with 

same mean and variance as 

Ji/n 1)/n I\z\<e z ^(ds, dz) and where Z\ = at + bB t + 

Jof\z\>e z N(ds,dz). 

2.3. Main result. We may now state our main result. 

Theorem 1. Assume that a : R i— ► M is bounded and Lipschitz continuous. Let 
e G (0, 1) and neH. There is a coupling between a solution {X t )t>o to (0) and an 
approximated solution {X p ,e , t s)t>a as in Subsection \2.2\ such that for all T > 0, 



E 



sup|X Pn(t) -X"; e (t) | 2 

[0,T] ' 



< C T (n- 1 + nS e (u)) , 



where the constant Ct depends only on T, <r, a, 6, 77i2(V). 

The first bound n" 1 is due to the time discretization (Euler scheme), and the second 
bound n8 e is due to the approximation of the increments of the Levy process. As 
noted by Jacod [5] , the first bound may be improved if there is no brownian motion 
6 = (but we have to work with some weaker norm) . 

We assume here that m^iy) < oo for simplicity: this allows us to work in L 2 . How- 
ever, we believe that Theorem [T] allows one to show that in the general case where 
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J m min(z 2 , l)v(dz) < oo, the family (n 1 + nS e (i/)) 1 sup[ 0T ] \X pn ^ — X™' e ^\ 2 is 

tight: decompose the Levy process Z t = Z\ + Zf, where Z 1 satisfies our assump- 
tions, and Z 2 is a compound Poisson process. Apply Theorem [1] between the jumps 
of Z 2 , and paste the pieces... this might be complicated to write, but the principle 
is very simple. 

2.4. Optimisation. Choose e = l/n. Then recalling that 8 € (v) < e 2 , we get 



E 



sup \X Pn(t) - X { ' t) | 

[0,T] F ' 



< Cr/n, 



for a mean cost to simulate (X™ '(t)")te[o,T] of order T(n + i 7 ' 1 /„(^)). 



• We always have F e {v) < Ce~ 2 , so that the cost is always smaller than Tn 2 . 

-o 



• If v(dz) ~ \z\ a dz for some a € (0,2), F e (i/) ~ e so that the cost is of 
order Tn" 1 ^ 1 '^. 

Still assume that v(dz) ~ \z\~ 1 ~ a dz, for some a € (0,2). When neglecting the 
small jumps, the mean cost to get a mean squared error of order l/n is of order 
rp n max{i,a/{2-a)) ^ gee [J2]) ; which is huge when a is close to 2. We observe that the 
present method is more precise as soon as a > 1. 

2.5. Discussion. The computational cost to get a given precision does not explode 
when the Levy measure becomes very singular near 0. The more v is singular at 
0, the more there are jumps greater than e, which costs many simulations. But 
the more it is singular, the more the jumps smaller than e are well-approximated 
by Gaussian random variables. These two phenomena are in competition, and we 
prove that the second one compensates (partly) the first one. 

Our result involves a suitable coupling between the solution (X t ) t > and its ap- 
proximation (X" ,e ) t >o. Of course, this is not very interesting in practise, since by 
definition, (X t )t>o is completely unknown. This is just an artificial way to estimate 
the rate of convergence in law, using a Wasserstein type distance. 

The simulation algorithm can easily be adapted to the case of dimension d > 2. We 
believe that the result still holds. However, the result of Rio [T^] is not known in 
the multidimensional setting (although it is believed to hold). We could use instead 
the results of Einmahl [5] . This would be much more technical, and would lead to 
a lower rate of convergence. 



3. Brownian approximation 

Consider the Levy process introduced in (JT|), consider x € K, a : K ^ K Lipschitz 
continuous, and the unique solution (X t ) t > to Recall consider a Brownian 
motion (W t )t>o, set 

(6) Z t = at + ^b 2 +m 2 {v)W u 

which has the same mean and variance as Zt. Let (Xt)t>o be the unique solution 
to 

(7) X t =x+ f a{X s _)dZ s . 

Jo 
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Theorem 2. Assume that a is Lipshitz continuous and bounded. Then it is possible 
to couple the solutions (X t )t>o to (0) and (X t )t>o to ^ in such a way that for all 
p > 4, allT > 0, all n > 1, recall (0j 



E 



sup \X t 

[0,T] 



XA 2 



< G 



,2/p-l 



m p {v 



,2/p . 



where Ct, p depends only on p, T, er, a, b, m 2 {y). 

If we only know that m^(y) < oo , then we choose n — [m^ (i^)~ 2 / 3 ], and we get, at 
least when m^v) < 1, E supr T i \X t — X t \ 2 < Cxm^u) 1 ' 3 . 

Consider a sequence of Levy processes (Z^) t >o with drift a, diffusion coefficient b 
and Levy measure v t , such that z 2 v e (dz) tends weakly to Sq. Then lim e ^o ^2(^e) — 
1, while in almost all cases, lim e ^o Tip(fe) = for some (are all) p > 2. 
Consider the solution to XI = x + a(Xl_)dZl. Then it is well-known and easy 
to show that (X^) t >o tends in law to the solution of a Brownian S.D.E. Theorem 
[2] allows one to obtain a rate of convergence (for some Wasserstein distance) . For 
example, we will immediately deduce the following corollary. 

Corollary 3. Assume that a is Lipshitz continuous and bounded. Assume that 
h '({\ z \ > e l) = f or some e € (0,1]. Then it is possible to couple the solutions 
{X t )t>Q to ([!)) and (X t )t>a to |7|) in such a way that for all r\ £ (0, 1), all T > 0, 



E 



sup \X t 

[0,T] 



XA 2 



< G 



.1-7) 



where Ct,tj depends only on rj, T, a, a, b, miiy). 

The original motivation of this work was to estimate the error when approximating 
the Boltzmann equation by the Landau equation. The Boltzmann equation is a 
P.D.E. that can be related to a Poisson-driven S.D.E. (see Tanaka [II]), while 
the Landau equation can be related to a Brownian S.D.E. (see Guerin [3J). In 
the grazing collision limit, the S.D.E. related to the Boltzmann equation has only 
very small jumps. However, many additionnal difficulties arise for those equations. 
Furthermore, we are able to prove our results only in dimension 1, while the kinetic 
Boltzmann and Landau equations involve 3-dimensional S.D.E.s 

4. Coupling results 

Consider two laws P, Q on R with finite variance. The Wasserstein distance W2 is 
defined by 

W 2 (P, Q) = inf {E [\X - Y\ 2 ] , C(X) = P, C(Y) - Q) . 

With an abuse of notation, we also write W 2 (X,Y) = W 2 {X,Q) = W 2 (P,Q) if 
C{X) = P and C{Y) = Q. We recall the following result of Rio [El Theorem 4.1]. 

Theorem 4. There is an universal constant C such that for any sequence of i.i.d. 
random variables {Yi)i>\ with mean and variance 9 2 , for any n > 1, 



c- 



n0 2 
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Here J\f(0, 6 2 ) is the Gaussian distribution with mean and variance 6 2 . Recall 
now ((3|). 

Corollary 5. Consider a pure jump centered Levy process (Y t )t>o with Levy mea- 
sure fi. In other words Yt = J J R zM(ds,dz), where M is a compensated Poisson 
measure with intensity dsfi(dz). There is an universal constant C such that 

Vi>0, Wl(Y t ,M(0,tm 2 (fi)))<C^\. 

m 2 (/zj 

Proof. Let t > 0. For n > 1, i > 1, write Y t n = n 1 ' 2 f{-^ )t/n / H zM(ds,dz), 
whence Y t = n -1 / 2 ^" The Y " are i i-d., centered, E[(F 1 ") 2 ] = tm 2 {n), and 

irE ( / / z 2 M (ds, dz) ] 



='/v 2 E 




z 2 M(ds,dz) + (£/n)m 2 (//) 



=n 2 [tm4(/i)/n + (tm 2 (/u)/n) 2 ] = ntm^^i) +t 2 m2(/i). 
Using Theorem |4j we get 



W|(y t ,^(0,tm 2 (M))) < C 
which concludes the proof. 



ntmi(fi) + £ 2 r7i 2 (/i) 
n£m 2 (/i) 



C 



m 4 (/x) 
m 2 (|U) ' 



□ 



This result is quite surprising at first glance: since the variances of the involved 
variables are im 2 (/i), it would be natural to get a bound that descreases to as 
t decreases to (and that explodes for large t). Of course, we deduce the bound 
Wf (Yt,jV(0, tm^in))) < Cmin(m4 (^i)/m 2 (/i), tm 2 (iu)), but this is now optimal, as 
shown in the following example. 

Example. Consider, for e > 0, /i e = (2e 2 ) _1 (<5 £ + S- e ), and the corresponding 
pure jump (centered) Levy process (Y" t e )t>o- It takes its values in eZ. Observe that 
TO 2(A i e) = 1 and m4(/i e ) = e 2 . There is c > such that for all t > 0, all e > 0, 
W 2 {Yt, 7V(0, t)) > cmin(t, e 2 ) = cmin(m4(/z e )/m 2 (/z e ), tm 2 (/J e )). Indeed, 

• if t < e 2 , then P(F/ = 0) > e -^«( R *) = e ~* /e2 > 1/e, from which the lowcrbound 
W$(Y f e ,N(0,t)) >ct = cmm(i, e 2 ) is easily deduced; 

• if t > e 2 , use that W$(Y t e ,J\f(0, t)) > E[min„ eZ \t x l 2 G - ne\ 2 } = ffi[min„ eZ |G - 
nei" 1 / 2 ) 2 ], where G is Gaussian with mean and variance 1. Tedious computa- 
tions show that there is c > such that for any a € (0, 1], E[min ng z \G — na\ 2 ] > 
(a/4) 2 P(G e U» eZ [(n + l/4)o, (n + 3/4)o]) > ca 2 . Hence W|(F t e ,^V(0,t)) > 
ct(et~ 1/2 ) 2 = ce 2 = cmin(t, e 2 ). 



5. Proof of Theorem Q] 

We recall elementary results about the Euler scheme for (J2]) in Subsection 15.11 We 
introduce our coupling in Subsection 15.21 which allows us to compare our scheme 
with the Euler scheme in Subsection l5.3l We conclude in Subsection l5.4l We assume 
in the whole section that a is bounded and Lipschitz continuous. 
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5.1. Euler scheme. We introduce the Euler scheme with step 1/n associated to 
©. Let 



A? 



(8) 

(9) = X^ 1)/n =Xf /n + a(xr /n )Ar +1 (i>0). 

The following result is classical. 

Proposition 6. Consider a Levy process (Z t ) t > as in fip. For (X t ) t > the solu- 
tion to (0j and for (X^ n )i>Q defined in 

E sup\X Pn(t) -Xl {t) \ 2 <C T /n, 

JO.T] 

where Ct depends only on T,a,b,m2{v), and a. 
We sketch a proof for the sake of completeness. 

Proof. Using the Doob and Cauchy-Scharz inequalities, we get, for < s < t < T, 

,t x 2 , „„ x 2 



(10) E 



sup\X u -X s \ 2 

[s,t] 



< CE 



a I \a(X u )\du) +sup(6/ a{X u )dB u 

[s,t] 



sup 

[a,t] 



a(X u _)zN(ds, du) 



<C T J (a 2 + b 2 +m 2 (iy))\\o-\\ 2 00 du<C T (t-s). 

Observe now that Xp n M = a; +/o'" a (^"„(s)-)^ s ' Setting A™ = sup[ t ] \X Pn ^ — 
X; n{s) \ 2 , we thus get A? = sup m \ j^ (s \n{X u ^) - a(X^ {u) _))dZ u \ 2 . Using the 
same arguments as in (fTO)) . then the Lipschitz property of a and (fTO]) , we get 

EL4™] < C T 



Pn(t) 



(a 2 + b 2 + m 2 {v))E[{a{X s ) - a(X;\ {s) )) 2 ]ds 



< c T / e[(x s - x Pn(s) y + (x Pn(s) - x; n(s) y]ds 



<C T / (\s-p n (s)\+E[A^])ds. 
Jo 

We conclude using that \s — p n (s)\ < 1/n and the Gronwall Lemma. 



□ 



5.2. Coupling. We now introduce a suitable coupling between the Euler scheme 
(see Subsection I5.1[) and our numerical scheme (see Subsection 12. 2|) . Recall §5§ . 



Lemma 7. Let neN and e > 0. It is possible to build two coupled families of i.i.d. 
random variables (A")^>i and (A"' e )i>i, distributed respectively as in fSJ) and 
in such a way that for each i > 1, 

E[(A»-AH 2 ]<CT 6 M, 

where C is an universal constant. Furthermore, for all e > 0, all n €E N, all i > \, 

E[A™] = E[A^ e ] = Var[A?] = Var[A^ e ] = b2+ ™^\ 



s 
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Proof. It of course suffices to build (A™, A™' e ), and then to take independent copies. 
Consider a Poisson measure N(ds,dz) with intensity measure dsv(dz)ln z \< e \ on 

[0,oo) x {\z\ < e}. Observe that J * fi z i <e zN(ds, dz) is a centered pure jump Levy 
process with Levy measure v e {dz) — \\ z \< t v{dz). Then we use Corollary [5] and 
enlarge the underlying probability space if necessary: there is a Gaussian ran- 
dom variable G™' e with mean and variance m2(v e )/n = mi^{v)/n such that 



E 



I Io /n L\< t zI *(ds, dz) - Gj' e | 2 < Cm 4 (^)/m 2 (^) = CS e (u). 



We consider a Brownian motion (B t ) t >o and a Poisson measure N with inten- 



sity measure dsv(dz)l{\ z \ >e y on [0. 



00 x 



{\z\ > e}, independent of the couple 



(G™' e , J J*M< e zN(ds, dz)) and we set 

. A™ := a/n + bB 1/n + tf /n Jj^ zN(ds, dz) + J^ /n f [z>e zN(ds, dz), 
. A^ £ := a/n + bB 1/n + G?' e + J^ /n Jj z|>e zN(ds, dz). 



while A™' e 



Then A™ has obviously the same law as Z\i n — Zq (see ([I) and $' 
has also the desired law (see (g}). Indeed, bB 1 / n + G"' c has a centered Gaussian 

law with variance b 2 /n + m,i^[v)/n = & 2 £ , and a/n + j^ z ^ >e zN(ds,dz) — 



a n,e + Jq^ 1 i| 2 |> e zN(ds, dz). This last integral can be represented as in Finally 
E[(A r 1 l - A™' e ) 2 ] < E [| J 1/n J {z{<e zN(ds, dz) ~ G"' e | 2 ] < CS e (v), and the mean and 
variance estimates are obvious. □ 



5.3. Estimates. We now compare our scheme with the Euler scheme. To this 
end, we introduce some notation. First, we consider the sequence (A™,A™' e )i>i 
introduced in Lemma [71 Then we consider PQ/„)i>o and (^/n)i>o defined in Q 
and ([5]). We introduce the filtration J 7 "' 6 = cr(A^, A^' e ,fc < i), and the processes 
(with V™' e = 0) 



Y i - A i/n 



X 



V-' 



i-l 
a \ -> 



a(X- /n )-a(X-; n )), 



Ml' 



k=0 



Lemma 8. There is a constant C, depending only on a, a, b, 1712(1') such that for 
all N > 1, 



E 



sup IK 

i=0,...,N 



n,ei2 



< Cn$ £ (i/)(1 + C/n) N (1 + iV 2 /n 2 ) 



Proof. We divide the proof into four steps. 

Step 1. We prove that for all i > 0, E [|K"' e | 2 ] < Cn8 e (u){l + C/nf. First, 



EIIF/^ 2 ] +E[(,(^JAf +1 -a(^)A^) 2 ] 



■2E 



(aPQJA^ - a(I^)A^) = E[|F 4 "' e | 2 ] + ^ + J, 



Now, using Lemma[7Jand that (A™ +1 , AK) is independent of J-™' e , we deduce that 

2a 

n 



^ B,e W/n)-^re)) 



n 



n,e |2i 



since a is Lipschitz continuous. Using now the Lipschitz continuity and the bound- 
edness of a, together with Lemma [7j and the independence of (A™ +1 , A™^) with 
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respect to T™' c , we get 

1^ < CE[\Yr\ 2 Wl\f] + CE[(A^\ - A? +1 ) 2 ] < -E[\Yr\ 2 ] + C5 e (v) 



Finally, we get 



E[|y^| 2 ] < (1 + C/n)E[\Yr\ 2 ] + C5 e (u). 



Since F n,£ = 0, this entails that E[\Y^ e \ 2 } < C6 e (v)[l + (1 + C/n) + ... + (1 + 
C/n) 1 - 1 ] < Cn8 e {v){\ + C/n) 1 . 

Step 2. We check that for N > 1, E[sup ^ |T/"' e | 2 ] < Cn5 e (2/)(1 + C/n) N N 2 /n 2 . 
It suffices to use the Lipschitz property of a, the Cauchy-Schwarz inequality, and 
then Step 1: 



E 



sup \Vi 

1....JV 



iV' 



JV-l 



^ E w 



< <4 E n\rr\ 2 } 



N-l 



< C^n5 t {v){l + C/n) 



N 



Step 3. We now verify that (M"' e )j> is a (jF™' e )i>o-martingale. We have M^_\ — 
M™' £ = a(X™ /n )[A? +1 - a/n] - a(X^p[A^ - a/n]. The step is finished, since the 
variables A™ +1 — a/n and A™^ — a/n are centered and independent of ^ r "' e . 

5*iep ^. Using the Doob inequality and then Steps 1 and 2, we get 



E 



sup \M, 

i=0....,N 



n,e\2 



<C sup E[|M™' £ | 2 ] 

i=0,...,JV 

<C sup E[|F"< e | 2 ]+C sup E[|U/ ; 

i=0,...,N i=0,...,N 

< CnS e (v)(l + C/n) N (l + N 2 /n 2 ). 



But now 



E 



sup \Y L 

i=0,...,N 



n,e\2 



< C*E 



sup \M, 

i=0,...,N 



n,e 1 2 



■CE 



sup | V i 

i=0,...,N 



n,e\2 



which allows us to conclude. 

Let us rewrite these estimates in terms of X n and X n ' e . 



□ 



Lemma 9. Consider the sequence (A™, A"' e )i>i introduced in Lemma^ and then 
(X^ n )i>o and (-^T^)i>o defined in and j3p. For all T > 0, 



E 



n,e |2 



< C T n(J 6 (i/), 



where Ct depends only on T,a 1 b,m,2{y),o '. 

Proof. With the previous notation, sup[ T j \X™ ^ — = sup i=0 [ nT j 

Thus using Lemma [HJ we get the bound CnS e (v)(l + C/n)^ nT \\ + [nT] 2 /n 2 ) < 
Cn5 t {v)e CT {\ + T 2 ), which ends the proof. □ 



\Y 
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5.4. Conclusion. We finally give the 

Proof of Theorem [7J Fix n G N and e > 0. Denote by Q{du, dv) the joint law 
of (A",A"' e ) built in Lemma and write Q(du,dv) = Qi(du)R(u, dv), where 
R(u, dv) is the law of A™' e conditionnally to A™ = u. 

Consider a Levy process {Z t )t>a as in ([T]), and (X t )t>o the corresponding solution to 
([2]). Set, for i > 0, A" = Zij n — Z(i_i)/„, and consider the Euler scheme (X^ n )i> 
as in ([5]). For each i > 1, let A"' e be distributed according to R(A™,dv), in such a 
way that (A™' e )i>i is an i.i.d. sequence. Finally, let (X™^)i>Q as in ([5]). 

By this way, the processes (X t )t>o, (X^ n )i>o and (-X^)i>o are coupled in such a 
way that we may apply Proposition [6] and Lemma [9l We get 



E 



sup \X pn(t) -X n '^ t) \ 2 

[0,T] v ' 



< 2E 



sup \X Pn(t) -X n {t) \ 

[0,T] 



< C T [n _1 +nS £ (u)}. 



This concludes the proof. 



□ 



6. Proofs of Theorem [2] and Corollary [3] 

We assume in the whole section that a is bounded and Lipschtiz continuous. We 
start with a technical lemma. 

Lemma 10. Let (X t )t>o an d (X t )t>o be solutions to (TJ[) and Then for p > 2, 
for all t > 0, all h G (0, 1], 



E 



sup \X t ~X t0 \P 

[ta-to + h] 



<C p (h p / 2 + hm p {v)), E 



sup \X t - X ta 

[t ,t +h] 



< C p hP/ 2 , 



where C p depends only on p, c, a, b, miiy). 

Proof. It clearly suffices to treat the case of (X t )t>o- Let thus p > 2. 

Using the Burholdcr-Davics-Gundy inequality and the boundedness of c, we get 



E 



sup \X t ~X t0 \P 
[to,t +h] 



< C„E 



to + h 



\aa(X s )\ds 



C P E 



t +h 



b 2 (T 2 {X s )ds 



< C p h p + CphP' 2 + Cpl 



p/2 



ta + h 



C p E 



t +h 



p/2' 



a 2 (X s )z 2 N(ds,dz) 



to Jl, 



p/2' 



z 2 N(ds, dz) 



< CphP' 2 + CpEp p h /2 ] 
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where Ut — /J / R z 2 N(ds,dz). It remains to check that for t > 0, E[[/f^ 2 ] < 
C p (t p / 2 + tm p (v)). But, with C p depending on m-iiy), 

E[Uf /2 } = [ ds [ u(dz)E[(U s + z 2 ) p/2 - U V J 2 \ 



<C p ds u(dz)E[z 2 Uf/ 2 - 1 + \z\ p ] <C p E[t/f /2_1 ]ds + C p m p {v)t 
Jo Jr, Jo 

<C P f EpP/^ds + CpieP/' 2 - 1 + m p {v))t, 
Jo 

for any e > 0. Hence E[Uf^ 2 } < Cp{t p/2 ~ 1 t+m p {v)t)e c v t t' L by the Gronwall Lemma. 
Choosing e = t, we conclude that E[Uf /2 } < C p (t p l 2 + m p (u)t). □ 

Proof of Theorem^ We fix n > 1, T > 0, and p > 4. 

Sfep i. Using Lemma [5] (see also Lemma [7]) we deduce that we may couple two 
i.i.d. families (A")^>i and (A")j>i, in such a way that: 

• (A")i>i has the same law as the increments (Zi/ n — ^(i-i)/ n )i>i of the Levy 
process (J]); 

• (A")i>i has the same law as the increments (Zi/ n — Zu_i\/ n )i>i of the Levy 
process ((6]); 

• for all i > 1, E[(A™ — A") 2 ] < Cm±{y) (we allow constants to depend on miiy)). 
Step 2. We then set X ™ = X ™ = i, and for i > 1, A?, = 1)/n +ff(XJj_ n/JA? 



and X", n = ^("_i)/„ + CT PQi-i)/ n )^i • Using exactly the same arguments as in 
Lemmas H] and El we deduce that E sup[ 0T ] \X^, t -, ~ ^"„(t)l 2 — Crnm^v), where 
Ct depends only on T, a, a, 6, m^ty). 

Step 3. But (X™ ( 4 pt>o is the Euler discretization of (O, while (X™ ^)t>o is the Eu- 
ler discretization of 0. Hence using Step 2, Proposition[6]and a suitable coupling as 
in the final proof of Theorem [TJ E supr T ] \X Pn ^ — X Pn ^\ 2 < Cri^/n + nm^u)). 



Step 4- We now prove that E 



sup [0iT] \X t - X pn(t) \ 2 < Cr.p^/P- 1 + m» 2 /P). 



We set T z = sup [l/n (l+1)/n] \X t - X Pn{t) \ = sup [i/n>(i+1)/n] \X t - X i/n \. By Lemma 
HOI E[Tf] < C p [{l/n)P/ 2 + m v (v)jn\. Thus, since p > 2, 

2/p 

E sup |X t -X Pri(t) | 2 <E sup T 2 <E sup rf 

[0,T] J l,...,[nT] l,...,[nT] 

< C T , P n 2/p [(1/nf/ 2 + m p (i/)/J , 

which ends the step. 

Step 5. Exactly as in Step 4, we get E supr 0T i \X t — X Pn ^\ 2 < CT, P n 2 / p ~ 1 . 

Step 6. Using Steps 3, 4 and 5, we deduce that with a suitable coupling, we have 
E[sup [0:T] \X t - X t \ 2 } < C T . v (n 2 lP- x + mp(yflP + n' 1 + nm 4 (v)). □ 

We conclude the paper with the 
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Proof of Corollary^ Since v{{\z\ > e}) = 0, wc deduce that m p {v) < rri2(v)e p 2 , 
for any p > 2. Applying Theorem[2]and choosing n — [e -p /( p-1 )] ) we get the bound 

C Tp (e {1 ~ 2/p){p/{p ^ 1)] + e (p- 2 H 2 /p) + e 2- P /(p-i)^ < C T p(e 1_1/(p_1) + e 2_4/p ). 

Hence for r) e (0,1), it is possible to get the bound Cr,jje 1_,? , choosing p large 
enough. □ 

Acknowledgement. I wish to thank Jean Jacod for fruitfull discussions. 

References 

[1] S. ASMUSSEN, J. ROSINSKI, Approximations of small jumps of Levy processes with a view 
towards simulation, J. Appl. Probab. Volume 38, Number 2 (2001), 482-493. 

[2] U. Einmahl, Extensions of results of Komlos, Major, and Tusnady to the multivariate case, 
J. Multivariate Anal. 28 (1989), no. 1, 20-68. 

[3] H. Guerin, Solving Landau equation for some soft potentials through a probabilistic approach, 
Ann. Appl. Probab. 13 (2003), no. 2, 515-539. 

[4] N. Ikeda, S. Watanabe, Stochastic differential equations and diffusion processes, Second edi- 
tion. North-Holland Mathematical Library, 24. North-Holland Publishing Co., Amsterdam; 
Kodansha, Ltd., Tokyo, 1989. 

[5] J. JACOD The Euler scheme for Levy driven stochastic differential equations: limit theorems, 
Ann. Probab. 32 (2004), no. 3A, 1830-1872. 

[6] J. Jacod, A. Jakubowski, J. Memin On asymptotic errors in discretization of processes, 
Ann. Probab. 31 (2003), no. 2, 592-608. 

[7] J. JACOD, T. Kurtz, S. Meleard, P. Protter The approximate Euler method for Levy 
driven stochastic differential equations, Ann. Inst. H. Poincare Probab. Statist. 41 (2005), 
523-558. 

[8] J. JACOD, P. Protter Asymptotic error distributions for the Euler method for stochastic 

differential equations, Ann. Probab. 26 (1998), no. 1, 267-307. 
[9] J. JACOD, A.N. Shiryaev Limit theorems for stochastic processes, second edition, Springer- 

Verlag, Berlin, 2003. 

[10] J. KOMLOS. P. MAJOR, G. Tusnady, An approximation of partial sums of independent rvs 

and the sample df I., Z. Wahrsch. verw. Gebiete 32 111-131 (1975). 
[11] P. Protter, D. Talay The Euler scheme for Levy driven stochastic differential equations, 

Ann. Probab. 25 (1997), no. 1, 393-423. 
[12] E. RlO, Upper bounds for minimal distances in the central limit theorem, to appear in Ann. 

Inst. Henri Poincare Probab. Stat., 2009. 
[13] S. Rubenthaler, Numerical simulation of the solution of a stochastic differential equation 

driven by a Levy process, Stochastic Process. Appl. 103 (2003), no. 2, 311-349. 
[14] H. Tanaka, Probabilistic treatment of the Boltzmann equation of Maxwellian molecules, Z. 

Wahrsch. Verw. Gebiete 46 (1978/79), no. 1, 67-105. 

Nicolas Fournier, LAMA UMR 8050, Faculte de Sciences et Technologies, Univer- 
site Paris Est, 61 avenue du General de Gaulle, 94010 Creteil Cedex, France, email: 
nicolas . f ournier@univ-parisl2 . f r 



